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SYSTEM, DEVICE, AND METHOD 
FOR TIME-DOMAIN EQUALIZER TRAINING 
USING A TWO-PASS AUTO-REGRESSIVE MOVING AVERAGE MODEL 

5 

CROSS-REFERENCE TO RELATED APPLICATIONS 

The present application may be related to the following commonly-owned 
10 United States patent applications, which are hereby incorporated herein by reference 
in their entireties: 

United States Patent Application No. XX/XXX,XXX entitled SYSTEM, 
DEVICE, AND METHOD FOR TIME-DOMAIN EQUALIZER TRAINING USING 
AN AUTO-REGRESSIVE MOVING AVERAGE MODEL, filed on even date 
15 herewith in the names of Qingli Liu and Aleksandar Purkovic (Attorney Docket No. 
2204/A07); 

United States Patent Application No. XX/XXX.XXX entitled SYSTEM, 
DEVICE, AND METHOD FOR TIME-DOMAIN EQUALIZER TRAINING USING 
INJECTED OUT-OF-BAND NOISE, filed on even date herewith in the names of 
20 Aleksandar Purkovic and Qingli Liu (Attorney Docket No. 2204/ A08); and 

United States Patent Application No. XX/XXX,XXX entitled SYSTEM, 
DEVICE, AND METHOD FOR DETERMINING THE SAMPLING TIME FOR A 
TIME-DOMAIN EQUALIZER, filed on even date herewith in the names of Qingli 
Liu and Aleksandar Purkovic (Attorney Docket No. 2204/A36). 

25 

FIELD OF THE INVENTION 

The present invention relates generally to communication systems, and more 
30 particularly to time-domain equalization of a communication channel. 

BACKGROUND OF THE INVENTION 

35 FIG. 1 shows a typical communication system 100 having a local/central unit 

102 in communication with a remote unit 106 over a communication medium 104. 
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Generally speaking, the communication medium 104 supports bi-directional 
communications between the local/central unit 102 and the remote unit 106. For 
convenience, communications from the local/central unit 102 to the remote unit 106 
are said to be "downstream" communications, while communications from the remote 

5 unit 106 to the local/central unit 102 are said to be "upstream" communications. 

Thus, the communication medium 104 typically supports a downstream channel over 
which the local/central unit 102 communicates to the remote unit 106 and an upstream 
channel over which the remote unit 106 communicates to the local/central unit 102. 
The upstream channel and the downstream channel may share the same physical 

10 communication link or occupy different physical communication links. When the 
upstream channel and the downstream channel share the same physical 
communication link, the upstream channel and the downstream channel may occupy 
the same frequency band (e.g., analog modem channels) or different, typically non- 
overlapping, frequency bands (e.g., ADSL or cable modem channels). The upstream 

15 and downstream channels may be symmetric or asymmetric. 

Within the communication system 100, it is common for the upstream and 
downstream communication channels to have dispersive characteristics. Specifically, 
each channel has a particular impulse response that disperses signals carried over the 
20 channel by extending the effects of each signal over a period of time. In many cases, 
the dispersive nature of the channel causes various distortions of the signals carried 
over the channel, such as Inter-Symbol Interference (ISI), Inter-Carrier Interference 
(ICI), and other distortions. 

25 FIG. 2 A shows a representation of an exemplary transmit signal 210 that is 

transmitted over a dispersive channel, for example, by the local/central unit 102. The 
exemplary transmit signal 210 includes two symbols, Si (211) and S2 (212), that are 
transmitted over the dispersive channel with no inter-symbol delay. 

30 FIG. 2B shows a representation of an exemplary receive signal 220 that is 

received over the dispersive channel, for example, by the remote unit 106, when the 
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symbols Si (21 1) and S 2 (212) are transmitted over the dispersive channel with no 
inter-symbol delay. As shown in FIG. 2B, the transmitted symbol Si (211) is 
dispersed by the dispersive channel such that the received symbol Ri (221) overlaps 
the beginning of the symbol S 2 (212). This causes ISI between the symbols Si (211) 
5 and S2 (212) and therefore corruption of the symbol S2 (212). 

One way to avoid or reduce ISI is to add a sufficient amount of inter-symbol 
delay to the transmitted symbols so that the received symbols do not overlap. 

10 FIG. 3 A shows a representation of an exemplary transmit signal 310 that is 

transmitted over a dispersive channel, for example, by the local/central unit 102. The 
exemplary transmit signal 310 includes two symbols, Si (31 1) and S2 (312), that are 
transmitted over the dispersive channel with inter-symbol delay. 

15 FIG. 3B shows a representation of an exemplary receive signal 320 that is 

received over the dispersive channel, for example, by the remote unit 106, when the 
symbols Si (31 1) and S2 (312) are transmitted over the dispersive channel with inter- 
symbol delay. As shown in FIG. 3B, the transmitted symbols Si (311) and S2 (312) 
are dispersed by the dispersive channel. However, because of the inter-symbol delay 

20 in the transmitted signal, the received symbols Ri (321) and R 2 (322) do not overlap. 
As a result, there is no ISI between the symbols Si (31 1) and S 2 (312). 

While the inter-symbol delay added to the transmitted signal eliminates (or at 
least reduces) ISI, there are detriments to employing such inter-symbol delay. For 
25 one, the inter-symbol delay reduces the efficiency of the transmitted signal in that 
fewer symbols (and therefore less data) are transmitted over a particular period of 
time. Also, the inter-symbol delay can cause cross-talk between channels carried over 
a common physical communication link. Thus, inter-symbol delay may be 
impractical for certain applications. 



30 
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Another way to avoid or reduce IS I is to "shorten" the impulse response of the 
channel. This is typically done using a time-domain equalizer (TEQ) at the receiving 
end of the communication channel. The TEQ is a short Finite Impulse Response 
(FIR) filter that is used to time-compress (shorten) the impulse response of the 
5 communication channel. In addition to shortening the impulse response of the 
channel, the TEQ also tends to "flatten" the channel and amplify noise. The 
effectiveness of the TEQ has a direct impact on overall performance, and therefore the 
TEQ design and the TEQ coefficients must be carefully determined. There are 
typically different design considerations for time-domain equalization of the upstream 
10 and downstream channels. Also, whether the implementation platform is memory or 
^ processing power limited (or both) plays an important role in the TEQ design. 

The following references are hereby incorporated herein by reference in their 
l~i entireties, and may be referenced throughout the specification using the corresponding 
f v 15 reference number. It should be noted that the reference numbers are not consecutive. 

y [1] John A. C. Bingham, ADSL, VDSL and Multicarrier Modulation, John 

h Wiley & Sons, 2000. 

O 20 [21 J. S. Chow, J. M. Cioffi, and J. A. C. Bingham "Equalizer Training 

Algorithms for Multicarrier Modulation Systems," ICC 1993, May 1993, pp. 761-765. 

[3] D. D. Falconer and F. R. Magee, Jr., "Adaptive Channel Memory 
Truncation for maximum Likelihood Sequence Estimator/' B.S.TJ. Nov. 1973, pp. 
25 1541-1562. 

[5] D. T. Lee, B. Friedlander, and M. Morf, "Recursive Ladder Algorithms 
for ARMA Modeling," IEEE Trans. Automat Contr., vol. AC-27, No. 4, Aug. 1982. 
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[6] P. J. W. Melsa, R, C. Younce, and C. E. Rohrs , "Impulse Response 
Shortening for Discrete Multitone Transceivers," IEEE Trans. Commun., vol. 44, No. 
12, pp. 1662-1672, Dec. 1996. 

5 [7] N. Al-Dhahir, A. EL Sayed, and J. M. Cioffi, "Stable pole-zero modeling of 

long FIR filters with application to the MMSE-DFE," IEEE Trans. Commun., vol. 45, 
No. 5, pp. 508-513, 1997. 

[8] N. Al-Dhahir, A. H. Sayed, and J. M. Cioffi, "A high-performance cost- 
10 effective pole-zero MMSE-DFE," Proc. Allerton Conf. Commun., Contr., Computing, 
Sept. 1993, pp. 1166-1175. 

[10] Steven M. Kay, Modern Spectral Estimation: Theory and Application, 
Prentice Hall, 1988. 

15 

[1 1] Peter E. Caines, Linear Stochastic Systems, John Wiley & Sons, 1988. 

[13] N. Al-Dhahir and J. M. Cioffi, "A low complexity pole-zero MMSE 
Equalizer for ML Receivers," Proc. Allerton Conf. Commun., Control, Comput, pp. 
20 623-632, 1994. 



SUMMARY OF THE INVENTION 

25 In accordance with an aspect of the invention, the coefficients for a time- 

domain equalizer are determined using a two-pass auto-regressive moving average 
model. A communication channel is first modeled using pi poles and qj zeros to 
form a first shortened channel impulse response having a first approximation Hi(z) = 
Bi(z)/1+Ai(z), wherein q% is greater than a predetermined cyclic prefix length. A 

30 time-mirrored image of the first shortened channel impulse response is then formed, 
and the resulting time-mirrored image of the first shortened channel impulse response 
is modeled using po poles and qi zeros to form a second shortened channel impulse 
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response having a second approximation H 2 (z) = B 2 (z)/1+A 2 (z), wherein q 2 is less 
than or equal to the predetermined cyclic prefix length. The time-domain equalizer 
coefficients are determined by combining Ai(z) and A 2 (l/z) with an appropriate 
amount of delay. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects and advantages of the invention will be 
appreciated more fully from the following further description thereof with reference 
to the accompanying drawings, wherein: 

FIG. 1 is a network diagram showing an exemplary communication system 
having a local/central unit in communication with a remote unit over a 
communication medium; 

FIG. 2A is a representation of an exemplary transmit signal having two 
symbols that are transmitted over a dispersive channel with no inter-symbol delay; 

FIG. 2B is a representation of an exemplary receive signal that is received 
over a dispersive channel when symbols are transmitted over the dispersive channel 
with no inter-symbol delay; 

FIG. 3A is a representation of an exemplary transmit signal having two 
symbols that are transmitted over a dispersive channel with inter-symbol delay; 

FIG. 3B is a representation of an exemplary receive signal that is received 
over a dispersive channel when symbols are transmitted over the dispersive channel 
with inter-symbol delay; 

FIG. 4 is a representation of an exemplary transmit signal having two symbols 
that are transmitted over a dispersive channel with cyclic prefixing; 

FIG. 5 is a representation of an exemplary receive signal that is received over 
a dispersive channel when symbols are transmitted over the dispersive channel with 
cyclic prefixing; 

FIG. 6 shows the impulse responses of simulated downstream and upstream 
G.lite ADSL channels; 
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FIG. 7 is a block diagram showing the relevant logic blocks of an exemplary 
ADSL terminal unit in accordance with an embodiment of the present invention; 

FIG. 8 is a block diagram showing the major logic blocks of an exemplary 
Discrete Multi-Tone (DMT) transceiver in accordance with an embodiment of the 
5 present invention; 

FIG. 9 is a representation of a time-domain equalized channel; 
FIG. 10 shows the shortened channel impulse response delay obtained from a 
reference model as known in the art; 

FIG. 1 1 shows the shortened channel impulse response delay obtained 
10 heuristically for a reference model as known in the art; 

FIG. 12 is a logic flow diagram showing exemplary logic for determining 
[V\ TEQ coefficients in accordance with an embodiment of the present invention; 
,y j FIG. 13 is a logic flow diagram showing exemplary logic for determining 

fu TEQ coefficients in accordance with an embodiment of the present invention; 

15 FIG. 14 is a logic flow diagram showing exemplary logic for determining 

f'= : TEQ coefficients in accordance with an embodiment of the present invention; 
r- FIG. 15 is a logic flow diagram showing exemplary logic for detecting 

J" instability in accordance with an embodiment of the present invention; 
■f- FIG. 16 shows the compressed (shortened) impulse response obtained by a 

pi 20 reference model for simulated downstream and upstream G.lite ADSL channels; 

FIG. 17 shows the compressed (shortened) impulse response obtained by an 
auto-regressive moving average model of the present invention; 

FIG. 18 shows a break-out of a channel impulse response into a head, a body, 
and a tail; 

25 FIG. 19 shows the signal to total distortion ratio obtained from a reference 

model and from an auto-regressive moving average model of the present invention; 
FIG. 20 demonstrates a two-pass shortening model of the present invention; 
FIG. 21 shows the signal to total distortion ratio obtained from a reference 
model and from a two-pass shortening model of the present invention; 
30 FIG. 22 shows the shortened channel impulse response delay obtained from a 

two-pass shortening model of the present invention; and 
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FIG. 23 is a logic flow diagram showing exemplary logic for determining 
TEQ coefficients using a two-pass shortening model of the present invention. 

5 

DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT 

An embodiment of the present invention applies the multichannel Levinson 
algorithm for auto-regressive moving average (ARMA) modeling of the channel 

10 impulse response. A two-pass shortening model is used for obtaining TEQ 

coefficients for a spread channel, such as the upstream channel of an Asymmetric 
Digital Subscriber Line (ADSL) system. In this two-pass shortening model, the 
communication channel is modeled in two passes, first for the original channel, and 
then to the time-mirrored image of the result from the first pass. The TEQ parameters 

15 are obtained by cascading the filter corresponding to the denominator polynomial of 
the first pass modeling and the filter corresponding to the mirrored image of the 
denominator polynomial of the second pass modeling. Delay is introduced so that the 
resulting shortened channel is closer to the optimal range. 

20 Various TEQ design considerations are discussed herein with reference to an 

exemplary communication system that is generally known as Asymmetric Digital 
Subscriber Line (ADSL). However, it should be noted that ADSL is simply an 
exemplary embodiment, and many of the TEQ design considerations are applicable in 
other networking environments. 

25 

An ADSL communication system provides high-speed digital communications 
between a telephone company central office and a subscriber site (e.g., a customer 
home or business) over existing twisted-pair copper telephone wires (i.e., the analog 
local loop, which is referred to in ADSL parlance as the "subscriber loop"). More 
30 specifically, the ADSL communication system provides high-speed digital 

communications between a central ADSL terminal unit (ATU-C) at the telephone 
company central office and a remote ADSL terminal unit (ATU-R) at the subscriber 
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site that communicate over the subscriber loop. The communication channel from the 
ATU-C to the ATU-R is considered to be the downstream channel, while the 
communication channel from the ATU-R to the ATU-C is considered to be the 
upstream channeL The upstream and downstream channels usually occupy non- 
5 overlapping frequency bands on the subscriber loop, with the upstream channel at a 
lower frequency band than the downstream channeL The downstream channel 
provides more communication bandwidth than the upstream channel, and therefore 
downstream communications are faster than upstream communications. 

Discrete Multitone (DMT) has been adopted as the modulation of choice in 
several ADSL standards. DMT was adopted first by ANSI (T1.413 standard) and 
then by the ITU-T (G. 992.1 and G. 992.2, formerly G.dmt and GJite, respectively) 
for high-speed data transmission over twisted-pair copper wire. In DMT-based 
transmission, the available signal bandwidth in either the upstream or the downstream 
direction is divided into multiple parallel orthogonal subchannels. The information 
bit stream is broken into multiple groups, with each group used to select a quadrature- 
amplitude-modulation (QAM) constellation point, which in turn modulates a 
subcarrier in the specific subchannel. One of the key design considerations for a 
DMT modem is the use of the Inverse Discrete Fourier Transform (IDFT) to modulate 
the DMT signal and the Discrete Fourier Transform (DFT) to demodulate the DMT 
signal. The use of DFT/IDFT requires that orthogonality between the subchannels be 
maintained at the receiver. The DFT and IDFT are effectively performed by using the 
Fast Fourier Transform (FFT) and Inverse FFT (IFFT), respectively. 

25 Efficient equalization of such subscriber loops can be carried out in the 

frequency domain by using a single-tap complex equalizer for each subchannel, but 
only if orthogonality between subchannels is maintained. This means there can be 
neither inter-carrier interference (ICI) nor inter-symbol interference (ISI) in the DMT 
signal. Also, in order to use DFT to perform demodulation, the received signal must 

30 appear periodic within a block of samples. The approach used in the DMT 
transceiver to accomplish this periodicity is to prefix a block of samples 
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corresponding to a DMT symbol, with the last v samples of the same block. The 
name used for these particular v samples is the cyclic prefix (CP). 

FIG. 4 shows a representation of an exemplary transmit signal 400 that is 
5 transmitted over an ADSL channel, for example, by the ATU-C over the downstream 
channel. The exemplary transmit signal includes two symbols, Si (410) and S 2 (420). 
Each transmit symbol consists of T=2N+v samples, where N equals the number of 
DMT subchannels. Each transmit symbol (410, 420) includes cyclic prefix samples 
(41 1 , 421) and data samples (412, 422), where the cyclic prefix samples (41 1 , 421) 
10 are the last v samples of the data samples (412, 422). 

ui While the CP can make the received signal appear periodic, it also reduces the 

H; information-carrying capacity of each symbol. Therefore, in order to maintain high 
f ■ j transmission efficiency, the CP length should be short relative to the DMT symbol 
HI 15 length. For this reason, the value vis fixed in various ADSL recommendations, 

typically to one-eighth the number of DMT subchannels (i.e., v = N/8). An ADSL 
j: implementation typically utilizes 256 DMT subchannels (G.dmt) or 128 DMT 
j2 subchannels (G.lite) on the downstream channel, and therefore vis equal to 32 or 16 
y samples, respectively, on the downstream channel. An ADSL implementation 
20 typically utilizes 32 DMT subchannels on the upstream channel, and therefore vis 

equal to 4 samples on the upstream channel. 

In order for the CP to make the received signal appear periodic, the amount of 
dispersion imparted to a transmitted symbol by the communication channel cannot 

25 extend past the CP of the following transmitted symbol. In other words, the impulse 
response of the channel cannot cause the length of the received symbol to be greater 
than 2N+v samples, where N equals the number of DMT subchannels and v equals 
the CP length. This is demonstrated by example in FIG. 5, which shows a 
representation of an exemplary receive signal 500 in which the received symbol Ri 

30 (510), corresponding to the transmitted symbol Si (410), does not extend past the CP 
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of the transmitted symbol S2 (420). This puts a significant constraint on impulse 
response of the communication channel. 

Unfortunately, most subscriber loops are relatively long and have impulse 
5 responses that are much longer than the CP length. This can be demonstrated in an 
exemplary G.lite system that is simulated using ANSI loop #7 with G.lite transceiver 
filters. FIG. 6(a) shows the impulse response of a simulated downstream G.lite 
channel. FIG. 6(b) shows the impulse response of a simulated upstream G.lite 
channel. The impulse response length of each channel (approximately 125 samples 
10 for the downstream channel and approximately 40 samples for the upstream channel) 
is significantly longer than the CP length (16 samples on the downstream channel and 
4 samples on the upstream channel). Thus, out of 2N+V received samples, more than 
v samples will be corrupted due to dispersion, resulting in both ISI and ICI in the 
received signal. 

15 

Solutions to this problem have been sought in both frequency and time 
domains. The most common approach is to use a TEQ to compress the impulse 
response of the channel ideally to be equal to or less than the CP length. Then, 
orthogonality between subcarriers is ensured, and a simple (single complex tap per 
20 subcarrier) frequency domain equalizer (FEQ) can be used after the DFT. 

In order for the TEQ to effectively compress the impulse response of the 
channel, the TEQ is "trained" to the impulse response of the channel. Specifically, the 
TEQ in the ATU-R is trained to the impulse response of the downstream channel, and 

25 the TEQ in the ATU-C is trained to the impulse response of the upstream channel. 
Training involves, among other things, determining a set of TEQ coefficients that 
provide the required amount of compression. Training is accomplished by 
transmitting a training signal over the channel and determining the TEQ coefficients 
based upon the received training signal. Thus, the ATU-R trains its TEQ based upon 

30 a training signal transmitted by the ATU-C over the downstream channel, and the 
ATU-C trains its TEQ based upon a training signal transmitted by the ATU-R over 
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the upstream channel. The training signal typically consists of repeated 2N samples 
of the pseudo-random number sequence, with no CP added. 

FIG. 7 shows the relevant logic blocks of an exemplary ATU 700. Among 
5 other things, the ATU 700 includes a data port 702 that is coupled to a DMT 

transceiver 704, and training logic 706 that is also coupled to the DMT transceiver 
704. The data port 702 provides an interface 701 to a computer system, such as a 
communication network at the telephone company central office (in the case of an 
ATU-C) or a personal computer at the subscriber site (in the case of an ATU-R), 
10 enabling the ATU 700 to send and receive information over the computer system. 
The DMT transceiver 704 provides an interface 705 to the subscriber loop, enabling 
the ATU 700 to send and receive information from another ATU over the subscriber 
loop. The training logic 706 includes various logic for training the components of the 
DMT transceiver 704, including, among other things, logic for sending a training 
15 signal over the subscriber loop and logic for determining TEQ coefficients based upon 
a received training signal. 

FIG. 8 shows the relevant logic blocks of an exemplary DMT transceiver 704 
in which only the main functions are included. The DMT transceiver 704 includes a 
20 transmitter 810 and a receiver 830 that are coupled to the subscriber loop 705 through 
a hybrid 820. 

At the transmitter 810, transmit data (TxD) from the data port 701 is first 
processed by the TX data processing block 811, which performs such functions as 

25 scrambling, forward error correction (FEC) encoding, and bit loading. The resulting 
bit stream is then processed by the constellation generation and scaling block 812, 
which divides the bit stream into small groups, selects a QAM constellation point for 
each subcarrier (with each constellation point representing a frequency domain 
sample), scales each frequency domain sample, and extends each block of N 

30 frequency domain samples (where N is the number of subcarriers used for 
transmission) to N2-2N samples in such a way that the samples are complex- 
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conjugate symmetric around the iV-th sample. The samples are then processed by the 
IDFT block 813, which performs an IDFT and forms a DMT block (symbol) from N 2 
output time-domain, real-valued samples. Each DMT symbol is processed by the CP 
block 814, which adds the CP (i.e., the last v samples of the DMT symbol) to the 
5 DMT symbol to form a prefixed DMT symbol The prefixed DMT symbol is 
processed by block 815, which converts the DMT symbol into a continuous-time 
signal for transmission over the subscriber loop 705 via the hybrid 820. 

At the receiver 830, the continuous-time signal received over the subscriber 
10 loop 705 via the hybrid 820 is first processed by block 83 1, which filters and samples 
the received signal. The resulting samples are then processed by the TEQ block 832, 
which performs time-domain equalization and CP removal so that a DMT time- 
domain symbol with N 2 samples is obtained. The resulting samples are then 
processed by the DFT block 833, which performs a DFT on the time-domain symbol 
15 in order to transform the time-domain symbol into N 2 frequency samples (where only 
the first N samples are taken). These frequency samples are processed by the block 
834, which performs gain scaling and FEQ on a per-subchannel basis, slices the 
sample in each subchannel to the nearest constellation point, and decodes the 
constellation points to recover the transmitted information bits. The resulting bit 
20 stream is processed by the RX data processing block 835, which performs such 
functions as bit unloading, forward error correction (FEC) decoding, and 
descrambling to produce receive data (RxD) that is forwarded to the data port 702. 

As discussed above, the primary function of the TEQ is to "shorten" the 
25 transmission channel of a DMT signal so that the combined (shortened) channel has 
energy confined to the CP window. In order for the TEQ to equalize the channel 
effectively, the maximum length of the shortened channel impulse response following 
time-domain equalization is (v-f 1) samples, where vis the CP length. Additionally, 
the channel flattening and noise amplification effects of the TEQ must be considered. 



30 
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The problem of determining the TEQ coefficients for channel impulse 
response shortening is discussed with reference to FIG. 9. In FIG. 9, xk represents the 
transmitted signal, hk represents the channel impulse response of the subscriber loop, 
rtk represents various noise components, y* represents the combined signal (channel 
5 impulse response and noise) presented to the TEQ, m represents the impulse response 
of the TEQ, r k represents the combined (shortened) impulse response of the channel, 
A represents the delay introduced by the transmission channel excluding the delay of 
the target channel, and bk represent the target impulse response (TER) of the target 
channel. It is assumed that the channel impulse response (hk) includes all transmit and 

10 receive transceiver filters in addition to the actual local loop channel. The noise term 
(/I*) includes crosstalk, thermal, and quantization noise, as well as radio frequency 
interference (RFI). TEQ coefficients (w k ) are chosen such that the impulse response 
of the combined channel and TEQ filters (upper branch) matches as close as possible 
the lower branch. The TIR (b k ) has maximum length of ( v+ 1) samples. A known 

15 pseudo-random sequence xk with the period equal to the one DMT symbol interval is 
used for the TEQ training. 

The error sequence e/c can be expressed as: 

20 e k = w k * y k -b k *x fc _ A = w T y k - b T x k (1) 

where "*" denotes linear convolution and the bold characters represent vectors (as 
well as matrices) as follows: 

25 w = 0<9, w ],..., WM-if, where M is the number of TEQ coefficients, 
b = [bo, bi,..., b v ] T , where vis the CP length, 

7 

= [Xk-4 Xk-A-h'.*, Xk-d-v\ , 
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The TEQ tap coefficients can be obtained, for example, by finding the values 
for TEQ (w), TIR (b), and delay (A) that minimize the Mean Squared Error (MSE), 
which is defined as: 

5 MSE = E{e\ } = E{(w T y k - b r xj 2 ; = w r R^w + b r R Jb - 2w r R^b , 
where the correlation matrices are defined as: 



10 



25 



R„=E{x k x k J, R„ =E{y k y k } t R J yx = E{y k x k } = (R^ f (2) 



In order to avoid a trivial solution in which the TEQ and TIR taps are equal to 
zero, some constraints are typically imposed on TEQ or TIR coefficients. 

One commonly used constraint requires that the norm of either the TEQ taps 
15 or the TIR taps equals unity (i.e., w T w =1 or b T b = 1, respectively). In this case, the 
equalizer taps are obtained by solving the eigenvector problem [3]. 

Another commonly used constraint is the "unit tap constraint" (UTC), which 
requires one of the TIR taps to equal unity [13]. This can be written as b T e, = 1, 
20 where e z - is a unit vector having all elements zero except for the /-th element, which 
equals unity. For convenience, a model that uses such a constraint is referred to 
hereinafter as a "UTCb" model, and is used as a reference model for performance 
comparisons against results obtained from various embodiments of the present 
invention. 



UTCb Reference Model 



30 



Briefly, the UTCb model obtains the preferred TEQ coefficients by 
determining TEQ coefficients for each Ae[0, N 2 -v-l ] and selecting as the preferred 
TEQ coefficients the TEQ coefficients providing the overall minimum MSE. The 
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UTCb model determines TEQ coefficients for a particular A by performing the 
following steps: 

Step 1 - Applying second order statistics (2) and the selected value A, calculate 
5 the following ( v+ 1 ) by ( v+l ) matrix: 

O a = R - R^ R^R^ 



Step 2 - Find minimum diagonal element (i 0 ptJopt position) of: 

10 Step 3 - Calculate TIR coefficients and corresponding MSE: 

CQW i 



Step 4 - Calculate corresponding TEQ coefficients: 



15 The preferred TEQ coefficients are those that give the overall minimum MSE, 

i.e., (Wopt)minMSE. 

The UTCb model performs well for time-domain equalization of ADSL 
channels. One problem with the UTCb model, however, is that implementation of the 
20 UTCb model is computationally intensive. This is due, in part, because iterations are 
performed for each Ae[0, N2-V-I]. 



Fewer iterations can be performed without loss of performance by limiting the 
search range to AefAj, N2- v-1], where Ai is the sample order number (k) of hk where 



1 13847-2204-A32 (12124MD) 
04/04/00 



-17- 

the block of vrnost powerful consecutive samples begins. FIG. 10 shows the value Aj 
for a simulated downstream G.lite channel. 

Furthermore, very little performance is lost if the UTCb model is applied for a 
5 single, heuristically found value of A = A 2 = A\ + M/2, where M is the TEQ filter 
length. As shown in FIG. 11, the value A 2 provides an acceptably low MSE, although 
it may not be the optimal value for A overall. 

Two-Channel Auto-Regressive (TC-AR) Model 

10 

In an embodiment of the present invention, the TEQ design is based upon an 
auto-regressive moving average (ARMA) model, which in turn is converted into a 
two-channel auto-regressive (TC-AR) model through embedding. The approach of 
embedding an ARMA model into a TC-AR model is discussed in [5]. Certain 
15 advantages of using embedding are discussed in [6]. 

In the ARMA model, the input and output signals of discrete multitones 
(DMT) are associated with the parameters of a zero polynomial and a pole 
polynomial, respectively, thereby resulting in a pole- zero model. The ARMA (pole- 
20 zero) model of the channel results in a non-Toeplitz autocorrelation matrix for which 
there are no computationally efficient solutions. 

However, embedding an ARMA model into a TC-AR model produces a model 
consisting of poles only. With embedding, the input and output signals form a 

25 column vector of dimension two, or two channels, and correlation parameters related 
to the predication vector of the TC-AR model are represented by 2x2 matrices. Thus, 
the TC-AR model of the channel results in a Toeplitz autocorrelation matrix, the 
elements of which are 2x2 matrices rather than scalars. Application of this approach 
to the shortening of long feed-forward and feed-back FIR filters of a channel 

30 equalizer is discussed in [7] and [8], 
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The parameters of the TC-AR model can be solved recursively using the well- 
known Levinson, Wiggins, Whittle, Robinson (LWWR) algorithm, which is referred 
to hereinafter as the multichannel Levinson algorithm. A full derivation of the 
multichannel Levinson algorithm is presented in [9], [10], and [11], 

Briefly, the ARMA (pole-zero) model assumes that the transfer function of the 
unknown system, H(zh h&s the form: 



H( z ) = ^ = ^- ~Iki = H(z) (3) 

X(z) X(z) 

jt=l 



where p is the number of poles, g the number of zeros, and X(z) and Y(z) are the input 
and output of the unknown system, respectively. An equivalent representation in the 
time domain is: 

15 e k =y k -y k (4) 
where y k is the output sample and y k is the modeled output sample. 

The output sample y k is represented by y k = . The ARMA (pole- 

20 zero) model is used to model the input-output relationship, which is characterized by 
the physical channel impulse response {h k } . Starting with a (p 9 q) pole-zero model, 
where p is the number of poles and q is the number of zeros of the model, the 
modeled output sample, y k , can be written as: 



25 % =-a x Vi a p%- P + b o x k+kXk-i + ■" + b q x k ^ (5) 
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For the cases where p > q, the objective is to estimate the ARMA 
parameters {a 1 ,- , s^»6o»* , *^ fl } based upon the knowledge of the second-order 
statistics of the input and output sequences { y k , x k } . Denoting these estimates at the 
jth recursion (1< j< p) by {a/ ? *-*,a ; ; ,^,--*fe/_^} where S = p-q, the output 
5 sample y k can be represented as: 

y* = h + 4 = y*-i - 4 y*- 2 a j y + & o + + - ■ - + b U x ^j-s + 4 > ( 6 ) 

where e } k is the jth-order residual error sequence and b( = 0 for i < 0. At the yth 
;;So recursion, an estimate of the ARMA model is determined in which the difference 
n between the number of poles and zeros is also 8 , i.e., equal to the desired difference 

u p-q- 

U Estimating the parameters of the ARMA model described in (5) needs the 

j 15 correlation matrix. In this case, it is a non-Toeptitz matrix for which there exists no 
|: computationally efficient solution. Instead, the method of emdedding an ARMA 

model into a two-channel AR model is used to generate a model consisting of poles 
2 only [5], [7], [8]. The correlation matrix then takes a Toeplitz form, and the 

multichannel Levinson algorithm can be used to solve the model parameters. 



20 



25 



Assuming that {x k } is white (which is approximately true for the time-domain 
DMT signal) and that y k _ { = y k ^ for 1 < i < p, (i.e., that the previous estimates are 



accurate), and defining a vector z k - 



„ x k+8 m 

converted to the two-channel-AR model as: 



, then the ARMA model of (6) can be 



=^!vi sJaw -®fr*-s ©K-; + 



el 



= 2 t + e t , 
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in which the forward prediction coefficient matrices of the linear predictor i k are 
defined as: 



-0/ = 



0 
0 



0 
0 



l<i<8-\ 
8<Li<j. 



(7) 



Applying the orthogonality principle, which deals with the real signal only and 
states that E[(z k - z k )z T k _ i ] = 0 (where / = 1 ,2,. . .J and T denotes transpose), results 
in: 



10 E 



m-1 



(8) 



15 



or equivalently: 



R(i) + 0/R(/-1) +— +9jR(f- y') = 0, l<i< j 



in which the autocorrelation matrix is defined as: 



(9) 



R(0 = £k.zLl= 



R„{i) R yx (i-5) 



= R r H). 



(10) 



20 



Combining (8), (9), and (10) for the case i = 0 results in: 



E[(z k -z k )zl] = E 



= R(0) + £e j,R(-/n) = E[e k e T k }= . (1 la) 



That is, 
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R«» + e 1 ; R(-i)+ • • • + g;.r(-7) = s{ 



(lib) 



Combining (9) and (1 lb) results in the yth-order two-channel AR model in 
5 matrix form, as follows: 



[ej ••• 0/ i] 



R(0) - R(-(;-l)) R(-y) 

R(7-l) ••• R(0) R(-l) 
R(j) R(l) R(0) 



=[o-o4 



(12) 



10 



Similarly, the equation for the backward prediction vector can be derived as: 



[i */ ... *{] 



' R(0) ••• R(-U-D) R(-j) 

RO-1) ■•• R(0) R(-l) 
R(j) - R(D R(0) 



= [z* o ••• o]. 



(13) 



Combining (12) and (13) results in: 



15 



I <Sf[ 



0/ I 



' R(0) ••• R(-(;-D) R(-;)' 

RO-1) R(0) R(-l) 

. R(j) •■• R(l) R(0) 



X» 0 ••- 0 



(14) 



In the next iteration (i.e., 1), (14) becomes: 



' R(0) ••• R(-j) R(-0' + D) 

RU) ••• R(0) R(-l) 
RO' + l) ••• R(l) R(0) 



0-0 



• (15) 
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The following are now defined: 







I 




'a -0 Lf" 


o ... 0 






I 


25 0 - 



(16) 



5 

and 



R(0) 



R(-;) R(-0" + D) n 



c /+] = 



Rtf) - R(0) R(-l) 
R(j + 1) •- R(l) R(0) 



10 



Combining (15) and (16) results in: 



0' +1 



e/ +1 i 



C ; ' +1 = 



(17) 



Rewriting the right side of (17) and also using (15) results in (18), as follows: 



15 



"K^[2* 0-0 /?] + [« 0-0 2^]' 
K>0-0I{] + [X5 0- 0£] 



k^ +1 [i */ •-■ <d; o]+[o g; ••• 0/ i] 
k* +3 [o e; •• 6/ i]+[i ••• *j 0] 



C ; + 1 



From (18), it can be seen that: 



20 



2* 0 
a 0 • 



■0/3 
• 0 2^ 



i ■■• 0 



0 0 



0/ 1 



(19) 



where: 
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a = 6 >R(1 ) + • • • + 0/R(;) + R( j + 1 ) , ( 20a ) 
y9 = OjR(-l)+--- + O/R(-;) + R(-0'+D). (20b) 

It can also be proven that: 

0 = 0?. ( 21 ) 

From (16), it can be seen that a + k£,,E* = 0 and fi + K b j+l Z f } = 0 , and 
therefore: 

Comparing the left-side of (17) with the right-side of (18), it can be seen that: 

0/' +1 = 0/ + K ' M <D \<i<) (23a) 
Of" 1 = + K* +1 0 1 < i < ;' ( 23b ) 

From (16), it can be seen that: 

Noting that a = - K' +1 X* and £ = - K* +1 X^ , (24a) and (24b) can be rewritten 



25 as: 



S/ +1 =(I-K;,K^, (25a) 
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E^^d-K^K^^- (25b) 

Similarly, by comparing the left-side of (17) with the right-side of (18), it can 
be seen that: 

e£}=K' (26a) 



~^j+l 

j+1 ~ j+1 



©E=K* (26b) 



Substituting (21) into (22b) and (26) into (22), (23), and (25) results in the 
10 equations of recursion (27), as follows: 



15 



20 



Recursions; 1 < j < p — l 

a = 0 jR(l ) 4- • • • + e/R( j) + RU + 1) (27a) 







(27b) 






(27c) 


0/ +1 =G/+0;::ioj_ ( . +1 


i<«<; 


(27d) 


o/ +1 =o/ + oj:i0;:. ;+1 


i<i<; 


(27e) 






(27f) 






(27g) 



25 



The initial conditions for the recursion (28) can be derived from the zeroth- 
order forward and backward predictors, as follows: 



Initial Conditions: 




0l=-R(l)(R(O)r 


(28a) 


Oj=-R(-1)(R(0))" 1 


(28b) 


i{= (i-e; o}) R (0) 


(28c) 


zf = a-*l ej)R(o) 


(28d) 



30 
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5 Upon the completion of recursions, the TEQ coefficients are obtained directly 

from the relationship shown in (7), that is, the first entry of matrix 0f , 1 <i< p , 



w n =a£ =1, 



w i =o/'=efai) 



i<i< p 



(29) 



=510 



For the case where p = q, that is, S = 0 , the TC-AR model takes a slightly 



different form. Defining z k = 



, the input-output relationship is now: 



z * 



i K 
o 1 



(30) 



115 where -0/ = 



-at 
0 



0 



The same recursions of (27) with the initial conditions of (28) can be used to 
estimate the model parameters. At the end of recursion, we have 
{af,---,a p p ,b* ,---,b p p } estimated with b£ determined as b p a =Z f p Q,2)/R„(Q). The 
20 TEQ coefficients are obtained in the same manner as in (29). 



For the case where p<q, the roles of the input and output need to be 
interchanged in order to use the recursions of (27) to estimate the model parameters. 
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Specifically, the two-channel vector z k = 
8 = # - p. The TC-AR model is then: 



is defined, where in this case 



A, , ft 



(31) 



where: 



-e/ = 



-&/ o 

0 0 

0 



I<z<£-1 



8<i< j. 



(32) 



Since is defined differently, the correlation matrix now takes the form of: 



R yx (i+S) R„(j) 



= R r H). 



(33) 



The recursions of (27) with the same initial conditions of (28) are directly 
applicable to this case for estimating {a|,---, a q p , , ••-,£*} with b%=l. The TEQ 
coefficients are then obtained as: 



w . = af =-0«,(U) 0</<p 



(34) 



Estimating the TC-AR model parameters requires the computation of 
correlation matrices defined either in (10) for p > q or (33) for p < q . These 
correlation quantities can be computed directly from the received time-domain DMT 
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signal and the known input signal during training. In order to suppress noise, 
correlation sequences may be computed repeatedly over many DMT symbols and 
then averaged. Alternatively, the impulse response of the physical channel can be 
estimated from the received DMT signal and the known input signal via FFT/IFFT, 
5 and the impulse response can be directly assigned to the correlation sequences. 

As described before, the time-domain DMT input sequence is approximately 

L 

white. Therefore, RJft^SjSt R y ,Q) = S x ^h m h m ^;R yx (l)=S x h t =R^( r l) 9 where 
S x is the average energy of the input signal. 

10 

ARMA model order selection is an important component in generating an 
accurate approximation of the channel. Subscriber loops in typical carrier service 
areas have different lengths and configurations (such as the existence or non-existence 
of bridged-taps). In order to more precisely model the impulse responses of 
15 subscriber loops, models with different pole-zero orders may be used. However, the 
subscriber loop characteristics are not known to the ATU before a connection is 
established, so the model order for a particular subscriber loop is not known a priori. 
Therefore, the model order must be determined dynamically. 

20 In an ARMA model, the model order is determined by the number of zeros 

and the number of poles. The number of zeros (q) is determined according to the CP 
length, such that q < (v+1). The number of poles (p) is constrained by the complexity 
of the system. If either is too large, then TEQ performance degrades. Therefore, they 
should not be made arbitrarily long. 

25 

Various tests have been proposed in the literature for determining p and q 
from the output sequence y% or from its autocorrelation function R yy [10]. One 
approach, discussed in [7] and [8], utilizes a "brute force" approach. In this approach, 
maximum allowable complexity is assumed, which in turn sets the upper bound for 
30 values of p and q. Next, an exhaustive search is performed over different choices of 
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p and q, within the allowable range, until an acceptable approximation is found based 
upon some predetermined statistic, which is typically the mean squared error (MSE), 
although other statistics such as signal to total distortion ratio (STDR) or shortened 
signal to noise ratio (SSNR) may alternatively be used. However, it is important to 
5 emphasize that increasing p and/or q could result in a worse approximation, 

depending on the channel impulse response This "brute force" approach is too 
time consuming and complex for certain implementations. 

An embodiment of the present invention reduces the number of searches 
10 needed to obtain the TEQ coefficients by only performing searches for (p - q) - x. 
Typically, no search is performed for p = 0. This is possible because the parameters 
hi of a low-order TC-AR model can be obtained from those of a high-order model if the 
Ul order differential between the pole and zero polynomials remains the same. 

1^15 Consider the forward prediction coefficient matrix defined in (7) for the fixed 

fli value of S, that is, the order differential between the pole polynomial and the zero 
polynomial. In each new recursion, one more prediction coefficient matrix is 
obtained while the matrices obtained in previous recursions are updated. Therefore, 

if"' ":: 

J:; as in the single-channel AR modeling case, the prediction coefficients of a low-order 
;r:20 model can be obtained from a high-order model. To provide reasonable performance, 

the lowest order may be limited to the one at which the predetermined statistic (e.g., 

MSE) is below a predefined level. 

In such an embodiment, searching starts from a high-order model with a fixed 
25 <X for example, the model with 32 poles for the case p > q or 32 zeros for the case 
p<q . Then correlation matrices are computed according to the high-order model 
and the multichannel Levinson algorithm is applied in order to determine TEQ 
coefficients for the particular value of S. At each recursion j\ { a/ +1 ,- } is stored, 

with al +1 =1 for p>qox { a£ +1 ,'- 9 aj*l_ s } for p<q. 



30 
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Multiple sets of TEQ coefficients are obtained after completing all recursions. 
For each set of TEQ coefficients, the received time-domain DMT sample sequence is 
equalized by a TEQ with the number of taps equal to the number of coefficients in 
that set of TEQ coefficients. The output sequence is shifted back by a fixed amount 
5 and then transformed to frequency domain via FFT. The input signal information is 
removed from the frequency samples, and then the IFFT is performed. The resulting 
time-domain sequence with a length of 2N samples yields the impulse response of the 
shortened channel. A sliding window with a length of (y + 1) is then used to find the 

( v +1) most powerful consecutive samples. A shortened signal to noise ratio (SSNR) 
10 is then computed. The SSNR is defined as the ratio of impulse response coefficient 
energy within a predetermined window to that outside the window, as follows: 

Ul a+v 

III 5 i=a 

Ol SSNR = 

t:: 4 a-1 2N-1 

i=0 i=a+v+l 

$1 where a is the number of the first sample of the ( v + 1) most powerful consecutive 
O samples. 

25 After a SSNR is obtained for each set of TEQ coefficients, the set of TEQ 

coefficients providing the highest SSNR is selected as the final TEQ coefficients. 

Thus, an embodiment of the present invention computes TEQ coefficients for 
each (p-q) = x(p*0) 9 finds the combined (shortened) channel for each such (p- q)~ 
30 x, computes the SSNR for each such (p-q)=x 9 and selects the TEQ coefficients 
providing the best SSNR. 

FIG. 12 is a logic flow diagram showing exemplary logic 1200 for 
determining TEQ coefficients using the multichannel Levinson algorithm for ARMA 
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modeling, wherein the coefficients for a low-order model are obtained from a high- 
order model. Beginning at step 1202, the logic computes TEQ coefficients for a new 
(previously untested) value of (p - q) = x, in step 1204. The logic preferably does not 
test the case for p = 0. The logic then finds the combined (shortened) channel 
resulting from the TEQ coefficients, in step 1206. The logic then computes the SSNR 
for the shortened (combined) channel, in step 1208. If the logic has completed all (p - 
q) combinations to be tested (YES in step 1210), then the logic selects the TEQ 
coefficients providing the best SSNR, in step 1212. If the logic has not completed all 
(p - q) combinations to be tested (NO in step 1210), then the logic recycles to step 
1204 to test (p - q) combination. The logic 1200 terminates in step 1299. 

When used specifically for determining TEQ coefficients for ADSL channels 
(as opposed to other types of communication channels), the number of searches can 
be reduced even further by searching only (p - q) = x for values of x > 0 (i.e., p > q). 
This is possible because simulation results have shown that, at least for ADSL, the 
TEQ coefficients for values of q > p are worse than the TEQ coefficients for values of 
p>q. 

FIG. 13 is a logic flow diagram showing exemplary logic 1300 for 
determining TEQ coefficients using the multichannel Levinson algorithm for ARMA 
modeling, wherein the coefficients for a low-order model are obtained from a high- 
order model and wherein only combinations of (p - q) for p > q are tested. Beginning 
at step 1302, the logic computes TEQ coefficients for a new, previously untested, 
value of (p - q) = x, where p > q, in step 1304. The logic preferably does not test the 
case for p = 0. The logic then finds the combined (shortened) channel resulting from 
the TEQ coefficients, in step 1306. The logic then computes the SSNR for the 
shortened (combined) channel, in step 1308. If the logic has completed all (p - q) 
combinations to be tested (YES in step 1310), then the logic selects the TEQ 
coefficients providing the best SSNR, in step 1312. If the logic has not completed all 
(p - q) combinations to be tested (NO in step 1310), then the logic recycles to step 
1304 to test (p - q) combination. The logic 1300 terminates in step 1399. 
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The number of searches can be further reduced by searching only for/? >/w 
where /w is a predetermined minimum value for p (e.g., p > 8). This is possible 
because, practically speaking, the model order will generally have at least the 
predetermined minimum number of poles, and therefore there is no need to search 
model orders having fewer poles. 

FIG. 14 is a logic flow diagram showing exemplary logic 1400 for 
determining TEQ coefficients using the multichannel Levinson algorithm for ARMA 
modeling, wherein the coefficients for a low-order model are obtained from a high- 
order model and wherein only combinations of (p - q) for/? > q and/? >p m in are tested. 
Beginning at step 1402, the logic computes TEQ coefficients for a new, previously 
untested, value of (p - q) = x, where p > q and p > pndn, in step 1404. The logic 
preferably does not test the case for p = 0. The logic then finds the combined 
(shortened) channel resulting from the TEQ coefficients, in step 1406. The logic then 
computes the SSNR for the shortened (combined) channel, in step 1408. If the logic 
has completed all (p - q) combinations to be tested (YES in step 1410), then the logic 
selects the TEQ coefficients providing the best SSNR, in step 1412. If the logic has 
not completed all (p - q) combinations to be tested (NO in step 1410), then the logic 
recycles to step 1404 to test (p - q) combination. The logic 1400 terminates in step 
1499. 

In theory, the TC-AR model is guaranteed to be stable, since the multichannel 
Levinson algorithm is used to estimate the model parameters [8], In a practical 
implementation, however, the TC-AR model, and, equivalently, the pole-zero ARMA 
model, can be unstable because of numerical errors. 

It can be observed from the recursions of (27) that the computations of 
a, ©£} and 3>£j , are error-prone. This is because, in computing a , successive 
matrix multiplications and accumulations are carried out. Small errors in the 
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parameters obtained from the previous recursion manifest themselves into significant 
errors in the new recursion. The computations of Q£{ and , on the other hand, 

involve matrix inversion. As the diagonal entries of Z£ and Z * become small in 

magnitude, inversions of these matrices yield large- value entries. Again, any small 
5 error in the original matrix results in a large error in the inverted one. 

Simulations have confirmed that, with very high signal-to-noise ratio (SNR), 
an unstable model can be obtained with a floating-point implementation. Thus, it is 
desirable to detect such instabilities during recursions. 



10 



Consider the case p>q. From (11a), it can be seen that: 



E{eV) E{*lx M ] 



(35) 



15 The upper diagonal entry of represents the mean square error of forward 

prediction. As the recursion value j increases, the mean square error decreases 
monotonically for a stable model. 



20 



Also, (27f) can be rewritten as: 
Zj +1 = (I~ej:io£)2;. (36) 



From (7), it can be observed that 0£J has non-zero entries on the first row. It 
can also be shown that O j*j has non-zero entries on the first column. The product of 

25 these two matrices has a non-zero entry only at the upper diagonal location, which is 
the only non-zero eigenvalue of the product matrix. For convenience, this eigenvalue 
is denoted as \ . Thus, the ARMA model is stable (i.e., the pole polynomial has all 
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its roots inside the unit circle in the z-plane or in other words, it has a minimum-phase 
response) if: 

0<A [ <1. (37) 

Equivalently, the matrix (I - ©£j<J>£j ) shall be positive definite to ensure the 

stability of the model. Also, its eigenvalues shall not be larger than 1 . Once this 
condition is met, it can be seen from (36) that the mean square error decreases 
monotonically. 

Thus, in an exemplary embodiment of the invention, the training logic uses 
(37) to check whether or not the model is stable. Specifically, the training logic 
monitors the value of ^ during each recursion to verify that the 0 < X x < 1 condition 
is met. If the model becomes unstable during a particular recursion (i.e., the 
0 < Aj < 1 condition is not met), that recursion is terminated, since any values that 
would be obtained thereafter would be unreliable. 

FIG. 15 is a logic flow diagram showing exemplary logic 1500 for detecting 
instability during an iteration of the multichannel Levinson algorithm for ARMA 
20 modeling. Beginning at step 1502, the logic performs a recursion of the multichannel 
Levinson algorithm for ARMA modeling, in step 1504. The logic then examines the 
value of Aj , in step 1506. If 0 < A 1 < 1 (YES in step 1508), then the logic recycles to 
step 1504 for a next recursion. Otherwise (NO in step 1508), then logic terminates 
the iteration. The logic 1500 terminates in step 1599. 

25 

An alternative is to add a small constant to the upper diagonal entry of 
2^ whenever its value reaches a low threshold. In this way, the stability of the model 

can be ensured for the current recursion. Similar correction can be done to . 
Although this correction is useful for obtaining a higher-order stable model, instability 
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of the model with a much higher-order cannot be completely avoided. Furthermore, 
with the correction, the mean square error may not decrease monotonically around the 
correction point. 

5 The multichannel Levinson algorithm for ARMA modeling can be expressed 

in a slightly different form. In this summary of the multichannel Levinson algorithm, 
it is assumed that the number of poles (p) is greater than the number of zeros (q), i.e. 
5>0. The algorithm takes a slightly different form for cases p=q and p<q, as will be 
apparent to a skilled artisan. However, since the best results are expected when p>q 

10 for time-domain equalization of ADSL channels, the algorithm is not described for 
the p<q cases for convenience. 



15 



The error sequence e* from (1) can be rewritten as: 

i=0 f=0 



(38) 



Introducing the following definitions: 



20 











, A,=< 




_ X k-S _ 









a ; 0 
0 0 

«i -Ks 

0 0 



where S=p-q, (38) can be rewritten as: 
= 5X u w with A 0 =1 = 



e, = u t - m 



>=0 



1 0' 
0 1 



l<i<S 
S<i<p 



(39) 



(40) 



25 It should be noted that (40) is similar to the AR (all pole) equation except that 

the variables are not scalers. Now, following the procedure in deriving the recursive 
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Levinson algorithm for the linear predictor, one can derive its multichannel (two- 
channel in this case) version. The following analagous variable applicable to the 
multichannel case are defined: 



10 



A/, i = 0,1,.. .,j , j-th recursion forward predictor coefficients 

j 

e J k = ^Ta/u^., A^'=I, ;-th recursion forward prediction error 



i=0 

i+1 



j-th recursion mean square forward prediction error 
j-th recursion backward predictor coefficients 



w / - X B /-^i u ^ 9 B o = 1 ' i~ thl recursion backward prediction error 



E J b = E{ w/ ( w J k ) T } , y-th recursion mean square backward prediction error 



R(0 =£K(u,f} = £< 



' R^ii) R yx {iS) 



= R r H) 



15 



20 



The multichannel Levinson algorithm can now be described as follows: 



Initial conditions: 



B°=I 



E° =R(0) 



E°=R(0) 



AIgorithm,7=0, l,...,p-l: 



25 Using A/ ,B/ for i = 0,1,..., j, and ,E/ calculate : 
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a' = ^AjRU-i + l), A'=I 

= fs B^ +1 R(-/), B^=I (41) 

5 

Calculate forward and backward prediction coefficients and MS errors: 

Af +1 =A/+KpB;:_ i+]> i = l,2 j 

10 Aj:l=K} 
B 0 W =I 

B/ +1 =B/+K^Aj. ;+1 , i = l,2,...,j 

Ef = (I-K}K^)E^ 
15 Ef 1 =(I-KiK})Ei 

Variables a 7 ,p ; , K } f , and , are only auxiliary variables. The last two 

variables are also known as the partial correlation coefficients. After p iterations, the 
model parameters are determined from the forward prediction coefficients (39). 

20 

Since the recursive procedures of the multichannel Levinson algorithm are 
carried out with 2x2 matrix operations, the multichannel Levinson algorithm is 
computationally efficient, particularly in comparison to the UTCb model. The 
multichannel Levinson algorithm is also robust with respect to numerical errors, and 
25 requires relatively low memory storage. However, use of the multichannel Levinson 
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algorithm to solve the parameters of the TC-AR model implies certain compromises, 
as discussed throughout the remainder of the specification. 

In an embodiment of the present invention, the parameters of the pole 
5 polynomial are directly applied to the TEQ, which forms a FIR filter. Specifically, 
the description of the ARMA model (3) can be rewritten as: 



15 



25 



H(z)=-^- =T^T«™ C 42 ) 



10 If the TEQ, W(z), is the denominator polynomial in (42), that is: 

W(z) = l +A(z) (43) 
then: 



H(z)W(z) = H(z)[l + A(z)l «^L[1 + A(z)l = B(z) (44) 

l + A(z) 



In other words, the equivalent (combined) channel made up of the original 
channel H(z) and the TEQ, W(z), in cascade, can be approximated by an FIR filter 
20 with the taps equal to the coefficients of the numerator polynomial, B(z). 

Consequently, the length of the combined channel approximately equals q, the 
number of zeros in the ARMA model. Therefore, if an accurate modeling is achieved 
with q < v+ 7, then ISI and ICI are avoided, and orthogonality between subchannels is 
achieved. 



It should be noted from (43) that the first TEQ coefficient is always one. This 
provides a non-trivial solution, but offers, for the same TEQ length, fewer degrees of 
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freedom compared to the UTCb model. For example, in the TC-AR model, there is 
no way to conveniently choose A, as there is in the UTCb model. Also, with the first 
TEQ tap always being unity, the resulting combined impulse response, tends to 
have the bulk of its energy located around the same position as the bulk of the energy 
5 for the original channel, This is not always optimal, especially when the original 
channel impulse response is spread. For the purpose of effective modeling, the first 
Ai samples corresponding to the "flat" part of precursor were not taken into account. 

Consider, for example, the upstream and downstream channels in an 
10 exemplary G.lite system that is simulated using ANSI loop #7 with G.lite transceiver 
O filters and white background noise at -lOOdBm/Hz. Even though the upstream and 
:;T downstream signals occupy different frequency bands, rather sharp filtering is still 
+; required in order to completely separate upstream and downstream signals. Because 
{i: the transceiver filters are considered to be part of the overall channel, the channel 
j|I 15 impulse response strongly depends upon their design. Since an upstream channel 
- requires a steeper filter slope than a downstream channel, the total channel impulse 

X: response is more spread for the upstream channel than for the downstream channel. 

Consequently, a longer TEQ is used for the upstream channel than for the downstream 
p channel (e.g., 32 taps for the downstream and 64 taps for the upstream). 

20 

When modeled using the reference UTCb model with the heuristically found 
delay A 2 > the TEQ provides effective equalization for both the upstream channel and 
the downstream channel. FIG. 16(a) shows the compressed (shortened) impulse 
response obtained from the UTCb model on the simulated downstream G.lite channel, 

25 with the compressed (shortened) impulse response shown as a solid line and the 
original channel impulse response shown as a dotted line. FIG. 16(b) shows the 
compressed (shortened) impulse response obtained from the UTCb model on the 
simulated upstream G.lite channel, with the compressed (shortened) impulse response 
shown as a solid line and the original channel impulse response shown as a dotted 

30 line. 
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When modeled using the TC-AR model, the TEQ provides effective 
equalization for the downstream channel but not for the upstream channel. FIG. 17(a) 
shows the compressed (shortened) impulse response obtained from the TC-AR model 
on the simulated downstream Giite channel, with the compressed (shortened) impulse 
5 response shown as a solid line and the original channel impulse response shown as a 
dotted line. FIG. 17(b) shows the compressed (shortened) impulse response obtained 
from the TC-AR model on the simulated upstream G.lite channel, with the 
compressed (shortened) impulse response shown as a solid line and the original 
channel impulse response shown as a dotted line. 

10 

In order to better measure and evaluate the achieved performance of the TC- 
AR model, a signal to total distortion ratio (STDR) is calculated. The total distortion 
is calculated as the sum of the combined ISI/ICI term and the noise terms, as 
measured at the slicer input. STDR is calculated as a function of the subcarrier order 
15 number, and is defined as follows: 

S 

STDR k = — , (45) 

where S R k is the receive signal discrete power spectrum, S l k is the discrete power 
20 spectrum of combined ISI/ICI due to insufficient channel shortening, and N k is the 
discrete power spectrum of the noise component. Index k represents the subcarrier 
order number. It is assumed that the samples of a non-ideally shortened channel 
impulse response do not span more than one DMT symbol interval (N2 samples). 
Those samples are denoted by pk } and are divided into three groups: 

25 

h = [ho, h h „. ? h H ~if, head, 

b = [bo> bj,..., b v ] T , body (most powerful v+i consecutive samples), 
t = 0>, t T -if, tail, 
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so that: 



P = \pa>pi,—> PNi-if = [ho> hi,..., h^i >bo, bi,... t b v t 0 , tj,..., t T .if (46) 

5 with the constraint, H + v + / + T = N 2 , as shown in FIG. 1 8. A rotated version of p, 
p\ is used as well, and its definition is: 

P' = \p'o>p'2,—,p'if3-j] T =[bo, bi,..., b v to, tj. h h 0 , hi,..., h H -if (47) 



10 Then, the components of (45) are computed as follows: 



N 2 -\ 



4 Pfi 



i=0 
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where, S x is the transmit (flat) power spectrum, calculated from the power spectrum 
20 density specification provided in the ADSL recommendations (-40dBm/Hz for the 
downstream and -38dBrn/Hz for the upstream); S NJt = S N is the discrete noise 

power spectrum, assumed flat for the purposes of this discussion; W(e J{0 ) is the 
Fourier transform of the discrete TEQ filter; 2L is the number of segments used for 
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numerical evaluation of the integral, that was originally in the expression for No = 
2N, where N is the number of subcarriers in the observed system, as stated before. 

FIG, 19(a) shows the STDR obtained on the simulated downstream G.lite 
5 channel, with the STDR obtained from the TC-AR model shown as a solid line and 
the STDR obtained from the UTCb model shown as a dotted line. FIG. 19(b) shows 
the STDR obtained on the simulated upstream G.lite channel, with the STDR obtained 
from the TC-AR model shown as a solid line and the STDR obtained from the UTCb 
model shown as a dotted line. 

10 

It is clear to see from FIG. 19(b) that straight application of generalized 
ARMA modeling does not produce satisfactory results in the upstream case. This 
may be explained as a consequence of several unfavorable facts. First, the number of 
zeros in the model is restricted to be less than or equal to the cyclic prefix. In the case 

15 of the upstream direction, v = 4, whereas in the case of the downstream, v= 16 

(G.lite) or v= 32 (G.dmt). Therefore, the upstream channel modeling is even more 
restricted, and it cannot be simply compensated by increasing the number of poles. 
Second, modeling the upstream channel with larger spread suffers even more from 
lacking the possibility to model delay, A, associated with the optimum MSE, as is 

20 done in the UTCb model. 

Two-Pass Shortening Model for Modeling the Upstream Channel 

Therefore, in an embodiment of the present invention, the upstream channel is 
25 modeled in two passes, first for the original channel, and then to the mirrored 
(reversed in time with appropriate delay) image of the result of the first pass. 
Specifically, the multichannel Levinson algorithm is first applied to the bulk of the 
original channel (the precursor samples are preferably discarded) with slightly relaxed 
requirements for shortening, such that the number of zeros in the first pass, qj, is 
30 greater than the CP length (i.e., qi > v). The multichannel Levinson algorithm is then 
applied to the mirrored (reversed in time with appropriate delay) image of the result of 
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the first pass. In this second pass, however, the number of zeros, q 2j is not permitted 
to exceed the CP length (i.e., q 2 < v). The TEQ parameters are obtained by cascading 
the filter corresponding to the denominator polynomial of the first pass modeling and 
the filter corresponding to the mirrored image of the denominator polynomial of the 
5 second pass modeling. The number of poles in the first and second pass, pi and p 2 , 
respectively, are chosen such that pi + p2+ 1 = M, where M is the TEQ filter length. 



An exemplary embodiment of the invention utilizes the heuristically 
determined values pj « 2M/3, qi ~2 v,p 2 - M/3, and q2 = v for the two-pass 
10 shortening model. 



The result of the two-pass shortening can be expressed in the following form: 



Usenet <z) = H (z) [1 + A (z)] [1 + 4 (z' 1 )] z <p ^ l) (48) 

15 

where H (z) is the original channel, [1 + A x (z)] is the denominator polynomial of the 
first pass modeling, and [1 + 4 {z~ l )] z" (P2+1) is the denominator polynomial of the 

second pass modeling, reversed in time and delayed. H shomned (z) can be evaluated 
based on the following expressions: 

20 

H (z) « 1 Bl< f! \ > first pass (49) 

i+4(z) 

H mirror ^ = ^ } ^N 2 ^ ^ fe) = jj + ^ ^ ^ tQ ^ 



H fror (z) ^ B 2 ( Z ) ? second (5l) 

25 

By combining (49), (50), and (51), expression (48) can be evaluated as: 
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H shortened (z) = B 2 (z- 1 ) z~ N > z^ l) (52) 

The TEQ coefficients are obtained by combining Aj(z) and Azillz) using the 
appropriate delay element, as described below. 

5 

The two-pass shortening model is demonstrated in FIG. 20. FIG. 20(a) shows 
the impulse response of the original channel. FIG. 20(b) shows the shortened impulse 
response after the first pass modeling. FIG. 20(c) shows the mirrored image of the 
shortened impulse response obtained from the first pass modeling. FIG. 20(d) shows 
10 the shortened impulse response after the second pass modeling. FIG. 20(e) shows the 
q final combined channel obtained from the two-pass shortening model. 

4- In order to better measure and evaluate the achieved performance of the two- 

ry pass shortening model, a signal to total distortion ratio (STDR) is calculated. The 
L!i5 total distortion is calculated as the sum of the combined ISI/ICI term and the noise 
» terms, as measured at the slicer input. FIG. 21 shows the STDR obtained on the 
"21 simulated upstream G.lite channel, with the STDR obtained from the two-pass 
O shortening model shown as a solid line and the STDR obtained from the UTCb model 
P shown as a dotted line. When compared to the STDR obtained from the single-pass 
^ = io TC-AR model, as shown in FIG. 19(b), it is clear that the performance of the two-pass 
shortening model is substantially better than the performance of the single-pass TC- 
AR model for time-domain equalization of the ADSL upstream channel. 

It should be noted that the two-pass shortening model introduces delay as a 
25 side-effect of modeling the mirrored image in the second pass modeling. Specifically, 
the flat precursor delay of the final combined channel, A, can be approximately 
determined with reference to FIG. 20(e), with the assumption that the shortened 
channel impulse response length equals the number of zeros (q) in the applied ARMA 
model. 



30 
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A s = A i + (3i-22)+0>2+l) (53) 

Since the number of zeros in the first pass modeling, q h is greater than the CP 
length and the number of zeros in the second pass modeling, q 2 , is less than or equal 
to the CP length, it follows that qi > q 2 , and consequently that A s > Aj. This has the 
desirable effect of shifting the shortened impulse response closer to the optimal range, 
as shown in FIG. 22. While both the UTCb model and the two-pass shortening model 
obtain satisfactory delay values, it should be noted that the UTCb model selects a 
value for A, whereas two-pass shortening model causes delay by the procedure. 

FIG. 23 is a logic flow diagram showing exemplary logic 2300 for 
determining TEQ coefficients for time-domain equalization of the upstream channel 
using the two-pass shortening model. Beginning at step 2302, the logic applies the 
ARMA model with p } poles and qj zeros to the bulk of the original channel impulse 
response, in step 2304, in order to form a first shortened channel having a first 
approximation Hi(z) = Bi(z)/1+Ai(z). The logic then forms a time-mirrored image of 
the first shortened channel, in step 2306, and applies the ARMA model with p 2 poles 
and q 2 zeros to the first shortened channel, in step 2308, in order to form a second 
shortened channel having a second approximation H 2 (z) = B 2 (z)/1+A 2 (z). The logic 
then combines Ai(z) and A 2 (l/z) with appropriate delay to obtain the TEQ 
coefficients, in step 2310. The logic 2300 terminates in step 2399. 

It is worth noting that the concept of modeling the reversed in time precursor 
of the impulse response was explored in [7] and [8]. Justification there was sought in 
the fact that the ARMA modeling achieves better results in modeling a decaying 
signal. While the multichannel Levinson algorithm guarantees a stable solution (all 
poles inside the unit circle)', assuming adequate processing precision, the reversing in 
time operation (replacement of z with 1/z) puts all poles outside the unit circle. In 
other words, the ARMA model obtained this way is unstable. Such instability is a 
problem for modeling the reversed in time precursor of the impulse response, as 
explored in [7] and [8], because the purpose there was to synthesize the channel. 
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Fortunately, such instability is not a problem for the two-pass shortening model of the 
present invention, since there is no attempt to synthesize the channel, but rather to 
determine the denominator polynomial coefficients that are used for setting the TEQ 
filter, which is always stable. 

5 

Alternative Embodiments 

In an exemplary embodiment of the present invention, the training logic 
including the two-pass shortening model is implemented predominantly as 
10 programmed logic that is stored in a computer readable medium and executed by a 
Digital Signal Processor (DSP) or other processor within the ATU. The programmed 
logic may be implemented in any of a number of programming languages. The 
yi programmed logic may be fixed in a tangible medium, such as a computer readable 
l|| medium (e.g., a diskette, CD-ROM, ROM, or fixed disk), or fixed in a computer data 
LHis signal embodied in a carrier wave that is transmittable to a computer system via a 
S1J model or other interface device, such as a communications adapter connected to a 
^ network over a medium. The medium may be a tangible medium (e.g., analog, 
:f : digital, or optical communication lines) or a medium implemented with wireless 
Jz techniques (e.g., microwave, infrared, or other transmission techniques). The 
';:J20 programmed logic embodies all or part of the functionality previously described 
herein with respect to the system. 

Alternative embodiments of the invention may be implemented using discrete 
components, integrated circuitry, programmable logic used in conjunction with a 
25 programmable logic device such as a Field Programmable Gate Array (FPGA) or 
microprocessor, or any other means including any combination thereof. 

The present invention may be embodied in other specific forms without 
departing from the essence or essential characteristics. The described embodiments 
30 are to be considered in all respects only as illustrative and not restrictive. 
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We claim: 

1 . A channel shortening method comprising: 

determining a first shortened channel impulse response for a communication 
5 channel using a first channel modeling scheme; 

forming a time-mirrored image of the first shortened channel impulse 
response; 

determining a second shortened channel impulse response for the time- 
mirrored image of the first shortened channel impulse response using a second 
10 channel modeling scheme; and 

combining the first shortened channel impulse response with an inverse of the 
O second shortened channel impulse response to obtain a third shortened channel 
y i impulse response. 

ilks 2. The channel shortening method of claim 1, wherein determining a first 
H - shortened channel impulse response for a communication channel using a first 

channel modeling scheme comprises applying an auto-regressive moving average 
22 model using pi poles and qi zeros to form the first shortened channel impulse 
'% response having a first approximation Hi(z) = Bi(z)/1+Ai(z), wherein g 2 is greater 
CI20 than a predetermined cyclic prefix length. 

3. The channel shortening method of claim 1, wherein determining a second 
shortened channel impulse response for the time-mirrored image of the first shortened 
channel impulse response using a second channel modeling scheme comprises 
25 applying an auto-regressive moving average model using p2 poles and qi zeros to 
form the second shortened channel impulse response having a second approximation 
H 2 (z) = B 2 (z)/1+A2(z), wherein q 2 is less than or equal to a predetermined cyclic 
prefix length. 

30 4. The channel shortening method of claim 1 , wherein the first shortened channel 
impulse response has a first approximation Hi(z) = Bi(z)/1+Ai(z) and the second 
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shortened channel impulse response has a second approximation HbCz) = 
B2(z)/1+A2(z), and wherein combining the first shortened channel impulse response 
with an inverse of the second shortened channel impulse response to obtain a third 
shortened channel impulse response comprises combining Ai(z) and A2(l/z) with 
5 appropriate delay. 

5. The channel shortening method of claim 1, wherein the first channel modeling 
scheme is based upon a pole- zero model having pj poles and the second channel 
modeling scheme is based upon a pole- zero model having poles, and wherein the 

0 number of taps in a channel shortening filter equals (pj + p 2 +1)* 

6. The method of claim 1, wherein the communication channel is an Asymmetric 
Digital Subscriber Line (ADSL) upstream channel. 
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7. An apparatus comprising a time-domain equalizer for equalizing a 

' communication channel and training logic for training the time-domain equalizer 
based upon a training signal received over the communication channel, wherein the 
training logic is operably coupled to determine a set of coefficients for the time- 
5 domain equalizer using a two-pass auto-regressive moving average model. 

8. The apparatus of claim 7, wherein the training logic comprises: 

first channel modeling logic operably coupled to determine a first shortened 
channel impulse response for the communication channel; 
10 inversion logic operably coupled to form a time-mirrored image of the first 

shortened channel impulse response; 
C ; : second channel modeling logic operably coupled to determine a second 

J : - shortened channel impulse response for the time-mirrored image of the first shortened 
i; ; channel impulse response; and 

HL 15 coefficient determination logic operably coupled to combine the first 

shortened channel impulse response with an inverse of the second shortened channel 
O impulse response to obtain a third shortened channel impulse response. 

J: 9. The apparatus of claim 8, wherein the first channel modeling logic is operably 
O 20 coupled to apply an auto-regressive moving average model using pj poles and gj 
zeros to form the first shortened channel impulse response having a first 
approximation Hi(z) = Bi(z)/1+Ai(z), wherein qj is greater than a predetermined 
cyclic prefix length. 

25 10. The apparatus of claim 8, wherein the second channel modeling logic is 
operably coupled to apply an auto-regressive moving average model using p2 poles 
and q 2 zeros to form the second shortened channel impulse response having a second 
approximation H 2 (z) = B 2 (z)/1+A 2 (z)> wherein q 2 is less than or equal to a 
predetermined cyclic prefix length. 



30 
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1 1 . The apparatus of claim 8, wherein the first shortened channel impulse 
response has a first approximation Hi(z) = Bi(z)/1+Ai(z) and the second shortened 
channel impulse response has a second approximation H 2 (z) = B2(z)/1+A 2 (z), and 
wherein the coefficient determination logic is operably coupled to combine Ai(z) and 

5 A2(l/z) with appropriate delay in order to obtain the third shortened channel impulse 
response. 

12. The apparatus of claim 8, wherein the first channel modeling logic is based 
upon a pole- zero model having pi poles and the second channel modeling logic is 

10 based upon a pole-zero model having p 2 poles, and wherein the number of taps the 
time- domain equalizer equals (pi +p2 +1). 



13. The apparatus of claim 7, wherein the communication channel is an 
Asymmetric Digital Subscriber Line (ADSL) upstream channel, and wherein the 
apparatus is a central ADSL terminal unit. 
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14. A program product for training a time-domain equalizer based upon a training 
signal received over a communication channel, the program product comprising: 

first channel modeling logic programmed to determine a first shortened 
channel impulse response for the communication channel; 
5 inversion logic programmed to form a time-mirrored image of the first 

shortened channel impulse response; 

second channel modeling logic programmed to determine a second shortened 
channel impulse response for the time-mirrored image of the first shortened channel 
impulse response; and 
10 coefficient determination logic programmed to combine the first shortened 

channel impulse response with an inverse of the second shortened channel impulse 
" ■ response to obtain a third shortened channel impulse response. 

jii 15. The program product of claim 15, wherein the first channel modeling logic is 
fU5 programmed to apply an auto-regressive moving average model using pi poles and qj 

zeros to form the first shortened channel impulse response having a first 
£j approximation Hi(z) = Bi(z)/1+Ai(z), wherein qj is greater than a predetermined 
tr cyclic prefix length. 

20 16. The program product of claim 14, wherein the second channel modeling logic 
is programmed to apply an auto -regressive moving average model using pi poles and 
q2 zeros to form the second shortened channel impulse response having a second 
approximation H 2 (z) = B 2 (z)/1+A 2 (z), wherein q 2 is less than or equal to a 
predetermined cyclic prefix length. 

25 

17. The program product of claim 14, wherein the first shortened channel impulse 
response has a first approximation Hi(z) = Bi(z)/1+Ai(z) and the second shortened 
channel impulse response has a second approximation H 2 (z) = B 2 (z)/l-i-A 2 (z), and 
wherein the coefficient determination logic is programmed to combine Ai(z) and 
30 A 2 (l/z) with appropriate delay in order to obtain the third shortened channel impulse 
response. 
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1 8. The program product of claim 14, wherein the first channel modeling logic is 
based upon a pole-zero model having p 2 poles and the second channel modeling logic 
is based upon a pole-zero model having pi poles, and wherein the number of taps the 
time-domain equalizer equals (pi + p2 +1). 

19. The program product of claim 14, wherein the communication channel is an 
Asymmetric Digital Subscriber Line (ADSL) upstream channel. 
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20. A communication system comprising a first communication device in 
communication with a second communication device over a communication channel, 
wherein the first communication device is operably coupled to transmit a training 
signal to the second communication device over the communication channel, and 
wherein the second communication device is operably coupled to determine a set of 
coefficients for a time-domain equalizer from a received training signal using a two- 
pass auto-regressive moving average model. 
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ABS TRACT OF THE DISCLOSURE 

5 A system, device, and method for time-domain equalizer training uses a two- 

pass auto-regressive moving average model. A communication channel is first 
modeled using p 2 poles and qi zeros to form a first shortened channel impulse 
response having a first approximation Hi(z) = Bi(z)/1+Ai(z), wherein q } is greater 
than a predetermined cyclic prefix length. A time-mirrored image of the first 
10 shortened channel impulse response is then formed, and the resulting time-mirrored 
image of the first shortened channel impulse response is modeled using p 2 poles and 
q2 zeros to form a second shortened channel impulse response having a second 
;U 5 approximation HfeCz) = B2(z)/1+Az(z), wherein qi is less than or equal to the 

i predetermined cyclic prefix length. The time-domain equalizer coefficients are 
f K J5 determined by combining Ai(z) and A2O/Z) with an appropriate amount of delay. 
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